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Abstract 

We discuss three Hamiltonians, each with a central-field part Hq and a PT-symmetric 
perturbation igz. When Hq is the isotropic Harmonic oscillator the spectrum is real 
for all g because H is isospectral to Hq+ g"^/2. When Hq is the Hydrogen atom then 
infinitely many eigenvalues are complex for all g. If the potential in Hq is linear in 
the radial variable r then the spectrum of H exhibits real eigenvalues for 0 < g < gc 
and a PT phase transition at gc- 
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1 Introduction 

It is known since long ago that some non-Hermitian operators may exhibit 
real eigenvalues BE]. This fact remained a somewhat exotic mathematical 
subject till Bender and Boettcher |3] suggested that those operators may ex¬ 
hibit unbroken parity-time (PT) symmetry. From then on the problem quickly 
developed into a prolihc held of research [1] (and references therein). 
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In a roughly general way we may say that most of the studied problems are 
represented by Hamiltonian operators of the form H = Hq + XH', where Hq is 
parity-invariant PHqP = Hq and H' is parity antisymmetric PH'P = —H', 
where P is the parity operator. U X = ig is imaginary (where g is obviously 
real) then H is PT symmetric: PTHPT = H, where T is the time-reversal 
operator |S]. 

In the beginning, most of the models studied were mainly one-dimensional 
[31IHI6U7] in which case Hq only exhibits parity symmetry and its eigenfunctions 
are even or odd: but later the researchers began to look 

for multidimensional examples |8HT5]. It was suggested that space-time (ST) 
symmetry could be a suitable generalization of the PT one [16]. In this case 
SHqS = Hq and SH'S = —H', where S' is a unitary operator such that 
Sl = S~^ = S. Most of the effort was devoted to hnd new multidimensional 
non-Hermitian Hamiltonians with real spectra. 

In the multidimensional case Hq and H may exhibit more complex symmetry 
that is conveniently described by means of group theory. In this way Fernandez 
and Garcia [TTlfTS] and Amore et al [T^[20] found that some ST-symmetric 
Hamiltonians exhibit broken ST symmetry for all values of g. The main con¬ 
jecture was that ST symmetry may be unbroken for some values of g provided 
that S is the only member of a class in the point group for Hq [20]. This 
appeared to be the case when S = P. 

The purpose of this paper is the discussion of three PT-symmetric Hamil¬ 
tonians for which Hq = /2 + V{r) exhibits central-held symmetry and 

H' = z. The resulting Hamiltonian H exhibits cylindrical symmetry and may 
be viewed as a kind of Stark effect with imaginary electric held. In Section [2] we 
outline the main ideas of PT symmetry as well as a simple argument based on 
perturbation theory [T7II20] . In Section [3] we briehy discuss the general case. In 
sections m 0 and ini we show that the models with V{r) = r^/2, V{r) = —1/r, 
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and V{r) = r, respectively, exhibit completely different spectra. Finally, in 
Section [7] we summarize the main results of the paper and draw conclusions. 


2 Parity-time symmetry 

Let A = PT = A~^ be the antiunitary operator given by the product of the 
parity P and time-reversal T operators [SIE|- The Hamiltonian operator H is 
said to be PT symmetric if 

AHA-^ = H. (1) 

If 

HtP = Ei), ( 2 ) 

then 

AHi) = AH A-^ Alp = HAip = AEip = E*Aip. (3) 

If Alp = aip, a being a complex number, then we say that PT symmetry is 
unbroken and E = E*. It is not difficult to prove that |a| = 1. Fernandez and 
Garcia [21] found a case in which Aip ^ aip and still E is real. They proposed 
the supposedly more general condition HAip = EAip for the occurrence of 
real spectrum; that is to say, when ip and Aip are two linearly independent 
eigenfunctions of H with the same eigenvalue E. This situation does not take 
place unless the spectrum of H is degenerate. However, further analysis reveals 
that both conditions are equivalent. In fact, if we choose ip = Ciip+C 2 Aip, where 
C 2 = aci and cl = ac 2 , then Ap = ap. It is worth adding that none of these 
conditions is of practical utility to predict whether H will have real eigenvalues 
or not because one commonly ignores the effect of A on the eigenvectors of H. 
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Most of the examples studied so far are of the form 


H = Ho + \H\ 


(4) 


where 


PHoP = Ho, THqT = Ho, PH'P = -H', TH'T = H' 


(5) 


and \ = ig, where g is real. Since TXH'T = X*TH'T = —XH' then AHA = H. 
Some useful information on the spectrum of H is given by the perturbation 
series 


e = J2 e^^^x\ 


( 6 ) 


j=o 


because if at least one coefficient of odd order is nonzero then E 


is expected to be complex for sufficiently small g. In such a case the PT- 
phase transition [T3] takes place at the trivial Hermitian limit = 0. If we 
write H{X)'ip^ = E^{X)'4)^ then PH{X)'iprn = PH{X)PP'ip^ = H{-X)Pxpm = 
Em{X)P'il!rn- If and P'lpm are linearly dependent, then Em{—X) = Em{X) 
and all the perturbation corrections of odd order vanish; otherwise P'lpm = '4^n, 
= En{X) and we cannot draw a conclusion so easily. The latter case 
may only take place when the spectrum of H is degenerate. In many cases it 
suffices to calculate the simplest, straightforward perturbation correction of 
hrst order EH [171120] . 

3 Stark effect 

Consider the Hamiltonian operator 
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(7) 


where V{r) is spherically symmetric (depends only on r). The eigenfunctions 
of Ho = H{X = 0) 




“ -^ul Vulm^ 


( 8 ) 


are also eigenfunctions of the angular momentum operators H and 

1 = 0,1,..., m = 0,±1,... ,±l. 


( 9 ) 


In spherical coordinates the eigenfunctions can be factored as 


( 10 ) 


where Rui{r) is the radial part, z/ = 0,1,..is the radial quantum number 
and Y{^{9, 0) are the spherical harmonics. Since the eigenvalues of Hq do not 
depend on m they are at least (21 + l)-fold degenerate. 

The perturbation H' = z breaks the degeneracy of the spectrum of Ho but 
the states with m > 0 remain two-fold degenerate because the eigenvalues of 
H do not depend on the sign of the magnetic quantum number m. 

Since 

HIL = (-i)V!,“h. (11) 


and PzP = —the matrix elements 


v'Vm 
I m 



( 12 ) 


are zero when / — /' is even. The perturbation corrections of hrst order to the 
energy are given by the eigenvalues of the matrix with elements z'^l'^. 
We will discuss three examples in the subsequent sections. 
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4 Isotropic harmonic oscillator 


When 

'^(f) = 


( 13 ) 


the Schrodinger equation is exactly solvable and the eigenfunctions and eigen¬ 
values are given by 


712 rL3 {x^y^z') (a))<y9^2( 2 /)) 

Ek=(k + ^ - -A^, k = ni + n2 + n^, 
ni,n2,n3 = 0, 1 ,..., 


(14) 


where y:>n{(l) is an eigenfunction of the one-dimensional harmonic oscillator 

ZJ _ I 1^2 

2 dq'^ ' 2 ^ • 


Since 


Mnin^nsix^y, Z) ='llJn^n2n3{-X, -y, -z)* = (-x) (“1/) + A*) 

= {-l)’^ijn,n2n3{x,y,z), (15) 

then the PT symmetry is unbroken for all g which accounts for the fact that 
the eigenvalues in equation ffT4)) are real for all g. 

Although in this case the approximate analysis based on perturbation theory 
may appear to be unnecessary we carry it out anyway merely for comparison 
purposes. To begin with, note that = {-lf'tl)^^}^^^^{x,y,z). 

The perturbation correction of hrst order to a given energy level is given 
by matrix elements of the form 

mi m 2 m 3 _ LtM 

^ni n2 ns \rni 112113 
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that vanish for all degenerate states because fc = ni +n 2 + ^3 = mi + m 2 + m 3 . 
Therefore, = 0 for all the states of the PT Stark effect in the isotropic 
harmonic oscillator. This result is consistent with the form of the exact eigen¬ 
values ffTT|) that depend on g"^. 

There is another way to prove that the PT symmetry for this problem remains 
unbroken for all values of g. The proof is based on the fact that H can be 
written in terms of a similarity transformation of Hq\ 

H = UHoU-^ + Y’ ^ (1'^) 

where Obviously, Hq and UHqU~^ are isospectral [ 22 ] . 


5 Hydrogen atom 


The unperturbed eigenvalues for the Coulomb interaction 

V{r) = (18) 

r 

are given by 

~ Yy ’ ^ ^ ^ + 1- (19) 

Therefore, there are pairs of degenerate states l^r which I — V = 

z/' — 1 / is odd and the corresponding matrix elements (I13 are nonzero. 

In such cases, which for real A give rise to what is commonly known as linear 
Stark effect [23112^ . the perturbation correction of hrst order is nonzero and 
the eigenvalues of H are complex for g ^ 0. 

The Schrodinger equation for this problem is separable in parabolic coordi¬ 
nates and the exact calculation of the perturbation corrections in terms of the 
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parabolic quantum numbers ni = 0, 1 , ..^2 = 0, 1 , ... and m = 0, ± 1 ,... is 
straightforward [25]. It is customary to write the perturbation series 

OO 

= ( 20 ) 

j=0 

in terms of the quantum numbers n = ni + n 2 + \m\ + 1 and q = ni—n 2 [25] . 
All the coefficients of odd order vanish when q = 0 but the states with q ^ 0 
are expected to be complex when g 0. 

The argument based on perturbation theory just outlined is sufficient to con¬ 
clude that this model exhibits complex eigenvalues when g ^ 0 and that the 
PT phase transition [T3] takes place at the trivial Hermitian limit g = 0. Nev¬ 
ertheless, we will show some numerical results to illustrate the point. Here we 
choose the most efficient method of Benassi and Grecchi [26] that is based on 
the separation of the Schrodinger equation in squared-parabolic coordinates. 
Since the details of this approach have been given elsewhere [261128] . here we 
just show the results. Figure [Tjshows the real and imaginary parts of the lowest 
eigenvalues. It is clear that the PT phase transition takes place at the trivial 
Hermitian limit as already argued above. 

The remarkable difference between the spectra of this problem and the previ¬ 
ous one can be traced back to the symmetry of Hq. The general central-held 
model is invariant under the group 0(3) while, on the other hand, the hy¬ 
drogen atom is invariant under the group 0(4) [22] . Such higher symmetry is 
due to the conservation of the Runge-Lenz vector in the latter model. Thus, 
the higher symmetry of Hq appears to be the reason why the PT symmetry 
is broken for all g in the perturbed hydrogen atom. While the /c-th harmonic- 
oscillator eigenvalue is *-^+^^^*-^+^) -fold degenerate, the n-th eigenvalue of 
the hydrogen atom is n^-fold degenerate. The greater degeneracy of the 
latter model allows the appearance of nonzero matrix elements and 

nonzero perturbation corrections of hrst order. 



6 Linear potential 


As a nontrivial example we consider the linear potential 
V (r) = r. 


( 21 ) 


In this case we cannot solve the eigenvalue equation for Hq exactly but we can 
nevertheless calculate the perturbation correction of first order to any energy 
level because it is determined by matrix elements of the form 


ly I m' 
I m 



( 22 ) 


which vanish for all sets of quantum numbers as argued in Section [31 Therefore, 
^^i\m\ ~ 0 there is a chance that PT symmetry may be unbroken for 
sufficiently small g. 

We can calculate approximate eigenvalues by means of diagonalization of a 
suitable matrix representation of the Hamiltonian. For simplicity, here we 
choose the nonorthogonal Slater-type basis set 

B = 0), n, /, |m| = 0,1,... (23) 


Present numerical results show that this problem exhibits the usual spectral 
pattern common to most PT-symmetric Hamiltonians studied by other au¬ 
thors; that is to say, unbroken PT symmetry for 0 < g < g^. For sufficiently 
small values of g the eigenvalues are real. As g increases two eigenvalues ap¬ 
proach each other, coalesce at an exceptional point [30fl33] gi > gc becoming 
a pair of complex conjugate numbers for g > gi- This behaviour is illustrated 
by figures [2], [3] and m for m = 0,1, 2, respectively. Those results were obtained 
by diagonalization of the matrix representation of the Hamiltonian operator 
in the Slater basis set (1231) with a = 2. The irregular lines reflect errors in the 
calculation of the eigenvalues originated in the quasi linear dependence of the 
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basis set. This shortcoming of the present approach becomes more noticeable 
as the number of radial basis functions increases. Although our numerical re¬ 
sults are not extremely accurate and are restricted to the lowest eigenvalues 
for the reason just indicated, they appear to suggest that the smallest excep¬ 
tional point Qc may be nonzero and that there is a PT phase transition at such 
point. We think that a more accurate calculation is not necessary to illustrate 
the difference between this model and the other two ones discussed above. 


7 Conclusions 


In this paper we have discussed three Hamiltonians given by three differ¬ 
ent central-held Hermitian parts and the same non-Hermitian PT-symmetric 
perturbation. Although at hrst sight they appear to be similar, they exhibit 
completely different spectra. In the case of the isotropic harmonic oscillator 
the PT symmetry is unbroken and the spectrum is real for all g. The reason is 
that H and Hq are related by the similarity transformation flT7|) . On the other 
hand, the PT symmetry is broken for all g in the case of the hydrogen atom. 
Quite in between the linear radial potential appears to exhibit unbroken PT 
symmetry for all 0 < g < gc and a phase transition at some gc that we were 
unable to determine. 

The remarkable difference among the spectra of such seemingly similar Hamil¬ 
tonians is due to the symmetry of Hq. As bl general rule the higher the sym¬ 
metry of Hq the more likely the occurrence of complex eigenvalues and the 
Hamiltonian for the hydrogen atom exhibits the greatest symmetry by far. 
We have already discussed the effect of symmetry in earlier papers [TTfED] 
but we have not seen such a remarkable difference in the behaviour of the 
non-Hermitian Hamiltonians. 

In closing we want to stress the fact that perturbation theory provides a useful 


10 


hint about the nature of the spectra of a given non-Hermitian Hamiltonian. If 
a perturbation correction of odd order (we typically look for the first one) is 
nonzero then we know that the spectrum is complex for all values of g (or at 
least for sufficiently small g). If all the available perturbation corrections of odd 
order are zero then there is a chance of finding real spectrum for some values of 
g. Obviously, this case should be investigated by more accurate calculations. 
As the symmetry of Hq increases, then also increases the dimension of its 
eigenspaces and, consequently, the dimension of the matrix representation of 
the perturbation in those eigenspaces. As a result it also increases the chance 
of nonzero perturbation corrections of first order. 

References 

[1] E. Caliceti, S. Graffi, and M. Maioli, Perturbation theory of odd anharmonic 
oscillators, Commun. Math. Phys. 75 (1980) 51-66. 

[2] G. Alvarez, Bender-Wu branch points in the cubic oscillator, J. Phys. A 28 
(1995) 4589-4598. 

[3] C. M. Bender and S. Boettcher, Real Spectra in Non-Hermitian Hamiltonians 
Having PT Symmetry, Phys. Rev. Lett. 80 (1998) 5243-5246. 

[4] C. M. Bender, Making sense of non-Hermitian Hamiltonians, Rep. Prog. Phys. 
70 (2007) 947-1018. 

[5] C. E. Porter, Eluctuations of quantal spectra, in: C. E. Porter (Ed.), Statistical 
theories of spectra: fluctuations, Vol. Academic Press Inc., New York and 
London, 1965. 

[6] P. M. Fernandez, R. Guardiola, and M. Znojil, Strong-coupling expansions for 
the PT-symmetric oscillators V(x) = a{ix) + b{ix)‘^ + c{ix)^, J. Phys. A 31 
(1998) 10105-10112. 

[7] F. M. Fernandez, R. Guardiola, J. Ros, and M. Znojil, A family of complex 
potentials with real spectrum, J. Phys. A 32 (1999) 3105-3116. 


11 



[8] C. M. Bender, G. V. Dunne, P. N. Meisinger, and M. Simsek, Quantum complex 
Henon-Heiles potentials, Phys. Lett. A 281 (2001) 311-316. 

[9] A. Nanayakkara and C. Abayaratne, Semiclassical quantization of complex 
Henon-Heiles systems, Phys. Lett. A 303 (2002) 243-248. 

[10] A. Nanayakkara, Real eigenspectra in non-Hermitian multidimensional 
Hamiltonians, Phys. Lett. A 304 (2002) 67-72. 

[11] A. Nanayakkara, Comparison of quantal and classical behavior of PT-symmetric 
systems at avoided crossings, Phys. Lett. A 334 (2005) 144-153. 

[12] H. Bila, M. Tater, and M. Znojil, Comment on: ’’Comparison of quantal and 
classical behavior of PT-symmetric systems at avoided crossings” [Phys. Lett. 
A 334 (2005) 144], Phys. Lett. A 351 (2006) 452-456. 

[13] Q-H Wang, Level crossings in complex two-dimensional potentials, Pramana J. 
Phys. 73 (2009) 315-322. 

[14] C. M. Bender and D. J. Weir, PT phase transition in multidimensional quantum 
systems, J. Phys. A 45 (2012) 425303. 

[15] C. R. Handy and D. Vrincenau, Orthogonal polynomial projection quantization: 
a new Hill determinant method, J. Phys. A 46 (2013) 135202. 

[16] S. Klaiman and L. S. Cederbaum, Non-Hermitian Hamiltonians with space-time 
symmetry, Phys. Rev. A 78 (2008) 062113. See also: Erratum: Non-Hermitian 
Hamiltonians with space-time symmetry [Phys. Rev. A 78, 062113 (2008)], 
Phys. Rev. A 89 (2014) 039908(E). 

[17] E. M. Eernandez and J. Garcia, Non-Hermitian Hamiltonians with unitary 
and antiunitary symmetries, Ann. Phys. 342 (2014) 195-204. arXiv: 1309.0808 
[quant-ph]. 

[18] E. M. Eernandez and J. Garcia, PT-symmetry broken by point-group symmetry, 
J. Math. Phys. 55 (2014) 042107. arXiv:1308.6179v2 [quant-ph]. 

[19] P. Amore, E. M. Eernandez, and J. Garcia, Non-Hermitian oscillators with Td 
symmetry, Ann. Phys. 353 (2014) 238-251. arXiv: 1409.2672 [quant-ph] 


12 


[20] P. Amore, F. M. Fernandez, and J. Garcia, Is space-time symmetry a suitable 
generalization of parity-time symmetry?, Ann. Phys. 350 (2014) 533-548. 
arXiv: 1405.5234 [quant-ph]. 

[21] F. M. Fernandez and J. Garcia, Critical parameters for non-hermitian 
Hamiltonians, Appl. Math. Commun. 247 (2014) 141-151. arXiv: 1305.5164 
[math-ph]. 

[22] F. M. Fernandez, Non-Hermitian Hamiltonians and Similarity Transformations, 
Int. J. Theor. Phys. (2015)(in the press) arXiv: 1502.02694 [quant-ph]. 

[23] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- And Two-Electron 
Atoms, (Springer, Berlin, Gotingen, Heidelberg, 1957). 

[24] L. D. Landau and E. M. Lifshitz, Quantum Mechanics. Non-relativistic Theory, 
(Pergamon, New York, 1958). 

[25] E. M. Eernandez, Introduction to Perturbation Theory in Quantum Mechanics, 
(CRC Press, Boca Raton, 2000). 

[26] L. Benassi and V. Grecchi, Resonances in the Stark effect and strongly 
asymptotic approximants, J. Phys. B 13 (1980) 911-930. 

[27] E. M. Eernandez and J. Garcia, Comment on ’’Stark effect in neutral hydrogen 
by direct integration of the Hamiltonian in parabolic coordinates”, Phys. Rev. 
A 91 (2015) 066501. On recent calculations of resonances for the Stark effect in 
hydrogen. arXiv: 1505.01327 [quant-ph] 

[28] E. M. Eernandez and J. Garcia, Highly accurate calculation of the resonances 
in the Stark effect in hydrogen, arXiv:1506.05797 [quant-ph]. 

[29] M. Lieber, 0(4) Symmetry of the Hydrogen Atom and the Lamb Shift, Phys. 
Rev. 174 (1968) 2037-2054. 

[30] W. D. Heiss and A. L. Sannino, Avoided level crossing and exceptional points, 
J. Phys. A 23 (1990) 1167-1178. 

[31] W. D. Heiss, Repulsion of resonance states and exceptional points, Phys. Rev. 
E 61 (2000) 929-932. 


13 


[32] W. D. Heiss and H. L. Harney, The chirality of exceptional points, Eur. Phys. 
J. D 17 (2001) 149-151. 

[33] W. D. Heiss, Exceptional points - their universal occurrence and their physical 
significance, Czech. J. Phys. 54 (2004) 1091-1099. 


14 





Fig. 1. Real and imaginary parts of the lowest eigenvalues of the PT-symmetric 
Stark effect in hydrogen 



Fig. 2. Lowest eigenvalues with m = 0 for the potential V{r,z) = r + igz 




Fig. 3. Lowest eigenvalues with m = 1 for the potential V (r, z) = r + igz with m = 1 




Fig. 4. Lowest eigenvalues with m = 2 for the potential R(r, z) = r + igz 
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